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ABSTRACT 

We discuss theoretical scenarios on crossover between nuclear matter (NM) and quark matter 
(QM). We classify various possibilities into three major scenarios according to the onset of diquark 
degrees of freedom that characterizes color-superconducting (CSC) states. In the conventional 
scenario NM occurs at the liquid-gas (or liquid-vacuum at zero temperature) phase transition and 
QM occurs next, after which CSC eventually appears. With the effect of strong correlation, the 
BEC-BCS scenario implies that CSC occurs next to NM and QM comes last in the BCS regime. 
We adopt the quarkyonic scenario in which NM, QM, and CSC are theoretically indistinguishable 
and thus these names refer to not distinct states but relevant descriptions of the same physical 
system. Based on this idea we propose a natural scheme to interpolate NM near normal nuclear 
density and CSC with vector coupling at high baryon density. We finally discuss the mass-radius 
relation of the neutron star and constraints on parameters in the proposed scheme. 

Subject headings: Neutron star, Equation of state, Nuclear matter, Quark matter. Color superconductor 


1. Introduction 

Quantum chromodynamics (QCD) is one fun¬ 
damental theory in the Standard Model in which 
non-Abelian gauge fields (gluons) and fermions 
(quarks) interact non-perturbatively and QCD- 
related phenomena are characterized by an energy 
scale; Aqcd — 0.2 GeV. The QCD-vacuum has 
a rich structure with a variety of condensates and 
it dynamically confines any colored excitation. At 
finite temperature of T ~ Aqcd and low baryon 
density of ^ ^qcd! ultra-relativistic nucleus- 
nucleus collision at Relativistic Heavy-Ion Collider 
called RHIC and Large Hadron Collider called 
LHC, supported by the lattice-QCD Monte-Carlo 
calculations, has revealed quantitative properties 
of a novel strongly-correlated state of QCD mat¬ 
ter, hadronic gas at T < Aqcd and the quark- 
gluon plasma (QGP) at T > Aqcd- In contrast. 


when T Aqcd and ub ^ ^qcd’ knowl¬ 
edge from QCD is severely limited because the 
first-principle numerical method based on the im¬ 
portance sampling fails at finite baryon density. 
Therefore, it is important to impose constraints 
on theoretical uncertainties from the experimen¬ 
tal side such as the neutron star observation and 
the beam-energy scan program at RHIC and fu¬ 
ture heavy-ion facilities. The purpose of this paper 
is to import the state-of-the-art idea of color de¬ 
confinement phenomenon at high T and high ne, 
especially a special feature of dense and large-A^c 
QCD, into the neutron star phenomenology, which 
in turn provides us with useful constraints. 

In some limiting cases at low and high baryon 
densities we have reasonable understanding of cold 
QCD matter. At ub = (1 ~ 2)no (where uq ~ 
0.16 fm~^ is the nuclear saturation density), on 
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the one hand, one can utilize empirical knowl¬ 
edge from nuclear physics and well-developed the¬ 
oretical methods to analyze nuclear matter (NM) 
properties (Weinberg 1990; Epelbaum et al. 2009). 
Recent developments include not only the chiral 
perturbation theory (Fiorilla et al. 2012; Krueger 
ct al. 2013) and the chiral effective model (Las- 
towiecki ct al. 2015; Holt ct al. 2014) but also 
the functional renormalization group (Drews & 
Wcisc 2015). On the other hand, at asymptoti¬ 
cally high density ub ^ lOOng (or even smaller 
densities hopefully), perturbative QCD (pQCD) 
calculations are validated for the bulk quantities 
of quark matter (QM) as seen from the conver¬ 
gence of three-loop perturbative results (Freedman 
& McLcrraii 1977; Kurkela et al. 2010; Fraga ct al. 
2015). Besides, pQCD works well to describe the 
color-superconducting (CSC) states. In particu¬ 
lar the ground state of QCD in the high-density 
limit has been identified as the color-flavor-locked 
(CFL) state (Alford ct al. 1999b) (see also a review 
by Fukushima & Hatsuda (2011) and references 
therein). 

The most problematic is the QCD matter study 
in the intermediate density region: 2no < ttb < 
lOOno or 1.1 GeV < /tb < 3 GeV in terms of 
the baryon chemical potential. Neutron stars are 
unique cosmic laboratories to access such domains 
experimentally (Buballa ct al. 2014). Useful ob¬ 
servations include the mass-radius {M-R) relation, 
the cooling curve (Shternin ct al. 2011; Page ct al. 
2004, 2006, 2011; Tsuruta et al. 2002; Blaschke 
ct al. 2000, 2004), the surface and toroidal mag¬ 
netic fields (Duncan & Thompson 1992; Cardall 
et al. 2001; Olausen & Kaspi 2014; Kitamoto ct al. 
2014), the gravitational waves from the merger 
of binary neutron stars (Abbott ct al. 2009; Ho- 
tokczaka ct al. 2011, 2013), and so on. In the 
present work we will specifically pay our attention 
to the mass-radius relation. 

The interesting point of the M-R relation is 
that it has a one-to-one correspondence to the the 
QGD equation of state (EoS) through the Tolman- 
Oppenheimer-Volkoff (TOV) equation. The over¬ 
all size of neutron star radii is determined by the 
EoS around riB = (1 ~ 2)no (Lattinicr & Prakash 
2001, 2007), while the typical neutron star masses 
M > Mq (where Mq denotes the solar mass) are 
largely correlated with the pressure at the central 
cores with ns = (2 ~ 10)no. According to the 


discoveries of PSR J1614-2230 (Demorest et al. 
2010) and PSR J0348-I-0432 (Antoniadis ct al. 
2013) there should be neutron stars whose mass 
exceeds 2Mq, indicating that the QCD EoS must 
be very stiff at ub > 2no as compared to what 
was naively considered. 

The requirement of stiffness challenges conven¬ 
tional EoS’s beyond the nuclear regime. In typ¬ 
ical hadronic models the strangeness appears at 
riB = (2 ~ 3)no (Glcndeimirig & Schaffncr-Biclich 
1999, 1998; Tsubakihara ct al. 2010), which sig¬ 
nificantly softens the EoS (i.e., hyperon puzzle). 
Seminal works (Nishizaki ct al. 2002; Vidana et al. 
2011; Weissenborn ct al. 2012; Yamamoto ct al. 
2014; Lonardoni ct al. 2015) introduced repulsive 
interactions of strangeness to circumvent the hy¬ 
peron puzzle, and the results imply that the re¬ 
pulsion not only in two-body force of baryons but 
also in three- and more-body forces are necessary 
to pass the 2Mq constraint. While many-body 
forces including hyperons have been poorly con¬ 
strained, new data on the nucleon-hyperon scat¬ 
tering and hypernuclei at J-PARC (Nagae 2010; 
Tamura 2010; Nakazawa & Takahashi 2010) as 
well as the lattice simulations at the physical 
quark masses (Inoue ct al. 2012; Beane et al. 2013) 
will provide important clues to resolve the hyperon 
puzzle. 

In general contributions from many-body force 
of baryons starts growing rapidly at tt-b ~ 2no 
(Hcbcler et al. 2010). From the microscopic point 
of view strong many-body correlations imply that 
baryons should exchange many mesons and quarks 
^ 2no. With more and more exchanged 
quarks, the identity of isolated baryon should be 
diminished, and quark degrees of freedom gradu¬ 
ally take over physical degrees of freedom. Eventu¬ 
ally, baryons with the radii (0.5 ^ 0.8) fm begin to 
overlap with each other at nn ^ (5 ~ 10)no, lead¬ 
ing to the percolated quark matter (Baym & Ghin 
1976). This physics picture contrasts with other 
crossover scenarios and fills smoothly in a gap be¬ 
tween nuclear and quark matter with exchanged 
mesons and quarks. The matter in the confine¬ 
ment/deconfinement crossover domain should in¬ 
herit properties from both nuclear and quark mat¬ 
ter. The idea of quarkyonic matter, that was first 
recognized by McLerran & Pisarski (2007), is an 
important step toward correct understanding of 
the crossover, as we will closely discuss later. We 
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will further push this idea forward concrete imple¬ 
mentation to construct an EoS. 

In this paper we delineate the crossover sce¬ 
nario bridged by quarkyonic matter. As we em¬ 
phasize later, the most essential point in this 
scenario is that there exists an overlapping re¬ 
gion that can be described in terms of strongly- 
interacting baryons or strongly-interacting quarks 
equally. For a concrete realization we utilize an 
Nambu-Jona-Lasinio (NJL) type model includ¬ 
ing vector and diquark interactions (for a review, 
Buballa (2005)) which are generally density de¬ 
pendent. Applying the quark-hadron duality pic¬ 
ture we determine the low-density behavior of 
parameters by fitting the Akmal-Pandaripande- 
Ravenhall (APR) nuclear EoS (Akinal ct al. 1998). 
The high-density behavior is constrained by con¬ 
fronting the EoS with the 2Mq constraint. In 
this way we cast the M-R relation onto the quark 
model parameters, and we analyze the transitional 
behavior with increasing density. 

One robust conclusion from our analysis is 
that the interactions should remain large even 
at riB lOriQ to account for the existence of 
massive neutron stars. One might consider that 
the asymptotic freedom in QCD would rather 
validate weakly interacting picture just in per¬ 
colated quark matter, since the typical distance 
among quarks is small. Such intuition, however, 
does not always work unless soft gluon contribu¬ 
tions are cutoff. For instance, in high-T QCD 
matter magnetic gluons remain unscreened and 
non-perturbative contributions survive even for 
asymptotically high temperature (Hietanon et al. 
2009). Also QCD in two-dimensions as an asymp¬ 
totic free theory (’t Hooft 1974) has dense quark 
matter in which non-perturbative gluons survive 
to asymptotically high density, since in spatially 
one-dimension screening effects are not enhanced 
at finite density (Schon & Thics 2000; Bringoltz 
2009; Kojo 2012). In dense QCD, while the hard- 
dense-loop (HDL) type calculations at weak cou¬ 
pling predict the electric screening mass of the 
QCD coupling constant times quark chemical po¬ 
tential, ^ gsfJ-qj actually the soft region in the 
gluon polarization function is highly dependent on 
the non-perturbative physics near the quark Fermi 
surface and can deviate from the HDL results at 
qualitative level (Rischkc et al. 2001; Huang & 
Shovkovy 2004b; Fukushima 2005; Kojo & Baym 


2014). Therefore the significance of the screening 
effects is still an open question until first-principle 
calculation determines the phase structure. In 
this situation it should be useful to construct an 
argument in which non-perturbative gluons sur¬ 
vive in percolated quark matter, as discussed in 
the quarkyonic scenario. In this paper we try 
to address this issue from the neutron star phe¬ 
nomenology, and use the 2 Mq constraint to dis¬ 
cuss why non-perturbative gluons should survive 

to ns ~ IOuq. 

The treatment of the crossover EoS proposed in 
the present work has some overlap with preceding 
works. The attempts to directly interpolate pQCD 
and nuclear EoS have been made by Freedman & 
McLcrran (1978); Fraga ct al. (2015). Also Alford 
ct al. (2005) discussed that an EoS for quark mat¬ 
ter with CSC can mimic the APR EoS. In these 
studies the descriptions of quark matter are based 
on the extrapolation of the pQCD picture. On 
the other hand, recent studies on the crossover 
scenario have used more model dependent but di¬ 
rect descriptions to conceptualize strongly corre¬ 
lated quark matter for ns ^ (3 ~ 5)no, then in¬ 
terpolated the resulting quark EoS and nuclear 
EoS at 2no (Masuda et al. 2013b,a, 2015; 

Alvarcz-Castillo ct al. 2014; Hell & Wcise 2014; 
Kojo ct al. 2015; Kojo 2015). These studies, how¬ 
ever, do not manifestly treat microscopic dynam¬ 
ics in the interpolated domain, so one cannot ad¬ 
dress, for instance, how the strangeness changes 
in the crossover domain. In the present work we 
will take a more direct and concrete path to the 
microscopic description with the interplay among 
quarks, diquarks, and baryons. We particularly 
put our emphasis on the role played by diquarks, 
with which there emerges a natural classification 
scheme of crossover scenarios as we see below. 

This paper is organized as follows. In Sec. 2 
we classify the crossover scenarios and discuss the 
underlying physics. In Sec. 3 we define our model 
and elucidate how to treat the model parameters. 
In Sec. 4 we construct an EoS and extract the 
density dependence of the model parameters. In 
Sec. 5 we discuss the resulting M-R relations and 
finally we make concluding remarks in Sec. 6. 
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2. Quark degrees of freedom beyond nu¬ 
clear matter 

There is no rigorous order parameter for quark 
deconfinement as soon as dynamical quarks are 
included in the theory, irrespective of various the¬ 
oretical efforts (Fukushima 2003). The Polyakov 
loop serves for an only approximate order param¬ 
eter at high T, while we have no clue about even 
an approximate order parameter if the baryon den¬ 
sity is high (see also an attempt by Dcxheimcr & 
Schramm (2010) to introduce a parameter analo¬ 
gous to the Polyakov loop). Here, we discuss the 
general features of deconfinement crossover not re¬ 
lying on any models. 

2.1. Three scenarios 

There is no established understanding on how 
quarks can become dominating beyond the normal 
nuclear density where nucleons are the most rele¬ 
vant degrees of freedom. Here we sort out various 
theoretical speculations on QM into three clearly 
distinct categories. We particularly pay special at¬ 
tention to the onset of diquark degrees of freedom 
introducing a diquark binding energy Bd together 
with a baryon binding energy Bb- Then, we can 
express the baryon mass as Mb = 3Mq — Bb and 
the diquark mass as = 2Mq — Bd^ where Mq 
represents the constituent quark mass (Fukushima 
et al. 2015). 

• Conventional scenario (NM < QM < CSC): 
In CSC at weak coupling QM is assumed 
to overcome NM beyond some baryon den¬ 
sity. Once we admit it, we can see that 
the Cooper instability suggests strong cor¬ 
relation of diquarks in momentum space. 
This, however, does not necessarily imply 
that those diquarks are localized in space 
or bound states. In fact we usually have 
Bd < 0 in the BCS-type calculation, leading 
to the following ordering: Mb/5 < Mq < 
Md/2. Therefore, there should be an on¬ 
set for QM next to the liquid-gas transition 
of NM. Such an onset is implemented in 
quark models or a first-order phase transi¬ 
tion between NM and QM around /ie 3Mq 
(with minor in-medium corrections) is often 
assumed to construct an EoS. At higher 
baryon density, eventually, diquarks form a 
condensate in the CSC states. 


• BEC-BCS scenario (NM < CSC < QM): In 
the intermediate region of baryon density the 
weak-coupling study may be drastically al¬ 
tered. If the attractive interaction is strong 
enough, the Cooper pair could become lo¬ 
calized in space and then CSC is identi¬ 
fied as a Bose-Einstein condensate (BEC) of 
bound diquarks (Abuki et al. 2002), which is 
called a (quasi-)molecule in condensed mat¬ 
ter physics systems controlled by a Eeshbach 
resonance (Ohashi & Griffin 2002). Then, 
we should anticipate a crossover transition 
between the BCS and the BEC regimes as 
the coupling strength changes (Nishida & 
Abuki 2005; Kitazawa et al. 2008; Sun et al. 
2007). In the BEC regime where Bd > 0, 
the mass hierarchy should be reorganized as: 
Mb/3 < Md/2 < Mq. This is actually a 
situation in some descriptions using quark- 
meson-diquark models in which diquarks are 
considered as physical degrees of freedom 
(Eukushima ct al. 2015). 

• Quarkyonic scenario (NM ~ QM ~ CSC): 
It is also a speculative but logical possibility 
that interacting nucleons are, in principle, 
indistinguishable from matter out of quarks. 
One heuristic measure to characterize quark 
deconfinement is the mobility of hadrons and 
quarks (Karsch & Satz 1980). Then, quarks 
can always hop from one nucleon to the other 
via meson exchange. In this sense the mo¬ 
bility of quarks is never vanishing as long 
as nucleons are interacting. One can up¬ 
grade this hand-waving argument to a more 
precise formulation by taking the large-Ac 
limit, which makes clear the difference be¬ 
tween the high-T low-/rB situation and the 
low-T high-/iB situation. In the former case 
mesons become non-interacting objects in 
the large-A^c limit, so that quark deconfine¬ 
ment is well-defined by percolation of over¬ 
lapping wave-functions of mesons. In the lat¬ 
ter case at high baryon density, on the other 
hand, the interaction among nucleons is of 
0{Nc) leading to a large pressure of 0{Nc) 
that is comparable to the pressure of quark 
matter (McLcrran & Pisarski 2007). This 
is a natural consequence from the fact that 
nucleon interactions are induced by quark 
exchange, and thus the pressure of nuclear 
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matter should be sensitive to quark degrees 
of freedom. This important observation on 
the pressure of 0{Nc) of both nuclear and 
quark matter opens a third scenario that 
there should be a regime of duality in which 
nuclear matter can be equivalently described 
in terms of quarks and diquarks (i.e., the 
McLerran-Pisarski conjecture as advocated 
by Fukushima & Sasaki (2013)). Such a spe¬ 
cial regime in the intermediate baryon den¬ 
sity is called the quarkyonic regime. 

In this work we adopt this last picture of the 
quarkyonic scenario. There are three reasons why 
we consider that the quarkyonic scenario should 
be the most realistic: 

(1) Center symmetry is more and more badly bro¬ 
ken with increasing baryon density, and so the 
deconfinement phenomenon at high baryon den¬ 
sity should be even broader crossover than that 
at high temperature. Thus, a first-order phase 
transition to quark matter is quite unlikely. In 
many models the chiral phase transition at low 
temperature could be of first order, but once a 
reasonable amount of vector coupling is included, 
it also becomes smooth crossover (Kitazawa et al. 
2002; Sasaki et al. 2007; Fukushima 2008; Bra- 
tovic et al. 2013). Also, the CSC phase can be 
smoothly connected to the confined phase; for ex¬ 
ample SU(N) lattice gauge theories coupled to 
fixed-length Higgs fields have a crossover region 
between confined and Higgs phase (Fradkin & 
Shcnker 1979). Therefore, it is the most con¬ 
ceivable that nucleons are gradually taken over by 
more fundamental degrees of freedom. 

(2) It is highly non-trivial how to reconcile con¬ 
finement and the BEC-BCS crossover. To form 
bound states of diquarks, the inter-quark inter¬ 
action should be strong, and then, the interac¬ 
tion among any colored objects including diquarks 
must be strong as well. Then, all the colored 
objects are expected to form color-singlet bound 
states. Theoretically speaking, however, what we 
can observe in principle is only the gauge invariant 
quantity and so there is no way to judge whether 
quarks and diquarks are confined in color-singlet 
bound states or not. This at the same time means 
that we cannot exclude a mixture of quarks and 
diquarks even in the confined phase where the in¬ 
teraction is strong. 


(3) The hadron resonance gas model is known 
to be successful to reproduce the thermodynamic 
properties of hadronic matter near or even slightly 
above Tc for small /ie (for a recent review, see 
Ding et al. (2015)). This is a clear example of the 
duality between hadrons and quarks in a transi¬ 
tional regime. From this point of view we can say 
that the quarkyonic regime could be a high-density 
counterpart of the so-called strongly-correlated 
QGP (i.e., sQGP) established in RHIG experi¬ 
ments. This analogy between sQGP and quarky¬ 
onic matter is addressed by Fukushima (2014). 

2.2. Quark-hadron continuity 

From the point of view of the symmetry¬ 
breaking pattern we can give some more solid ar¬ 
guments to justify the quarkyonic scenario which 
is quite consistent with the GSG theory. It has 
been a closely investigated conjecture called the 
quark-hadron continuity that the superfluid nu¬ 
clear matter is a dual state of the CFL phase 
(Schaefer & Wilczek 1999; Alford ct al. 1999a; 
Fukushima 2004; Hatsuda et al. 2006). The point 
is that the superfluid nuclear matter and the CFL 
phase break chiral symmetry and U(1)b symme¬ 
try. The theoretical description based on the sym¬ 
metries is elegant and convincing only in the chiral 
limit, but the strange quark mass comparable to 
Aqcd introduces subtlety. In particular the CFL 
phase may not continue to the normal nuclear 
density due to the Fermi surface mismatch. Then, 
possibilities include a CFL state entering the gap¬ 
less region which is known to suffer instabilities 
against spatial modulation (Fukushima et al. 2005; 
Huang & Shovkovy 2004a; Fukushima 2005). If 
the interaction is sufficiently strong, the instability 
might be avoided (Gubankova et al. 2006). 

In this work we postulate a large enough value 
of the diquark coupling H, so that we do not 
have to cope with instability problems. Then, 
there should be a phase transition to the two- 
flavor super-conducting (2SC) phase down from 
the CFL phase as lowering the baryon density. 
The pure 2SC phase has diquark condensates of 
ru-gd and rd-gu quark pairs in color space of red 
and green and in flavor space of up and down. 
Within the two-flavor sector these combinations 
in flavor space are SU(2) singlets, and so none 
of SU(2 )l X SU(2)r is broken. Besides, U(1)b 
is broken but a modified U(l)g is kept unbro- 


5 


ken, where U( 1)0 is generated by a mixed charge 
B — B — 2Qf,. In fact, the diquarks ru-gd and rd- 
gu have B = 1/3 + 1/3 = 2/3, the electric charge 
Qf, = 2/3 — 1/3 = 1/3 in unit of e, and therefore 
B = 0. Hence, no global symmetry is broken in 
the 2SC phase. 

It is also possible to have not the pure 2SC 
phase but the coexisting 2SC phase with non-zero 
chiral condensates. In this case chiral symmetry 
is spontaneously broken by the chiral condensates, 
and the diquark condensate breaks no new sym¬ 
metry. Thus, the coexisting 2SC phase has a com¬ 
pletely identical pattern of the symmetry breaking 
as the hadronic phase. In other words we have no 
way to exclude the 2SC-type diquark condensate 
even in NM in the confined phase. Considering the 
strange quark mass, such continuity between the 
coexisting 2SC phase and NM should be a more re¬ 
alistic candidate than the CFL-NM continuity. In 
conclusion, the quarkyonic scenario based on de¬ 
confinement is supported from the chiral symme¬ 
try point of view by the presence of the coexisting 
2SC phase with non-zero chiral condensates. 

Here, let us make two remarks to clarify pos¬ 
sibly confusing points. The first point is the or¬ 
der parameter and gauge symmetry. In the mean- 
field approximation there can be a phase transi¬ 
tion between normal QM and the 2SC phase and 
the associated order parameter is the diquark con¬ 
densate. So, one might think that there should 
be also a phase transition from NM to the 2SC 
phase even though the symmetry has no differ¬ 
ence. This is true only approximately as long 
as one fixes the gauge and relies on the mean- 
field approximation. However, generally beyond 
the mean-field level, there is no gauge-invariant 
order parameter for the 2SC phase while one 
could construct higher dimensional order param¬ 
eters for the CFL phase (Rajagopal & Wilczek 
2000; Fukushima & Hatsuda 2011; Fukushima 
2004). This “non-existence” of the 2SC order pa¬ 
rameter strongly suggests crossover between NM 
and the 2SC phase, which is reminiscent of the 
non-existence of the exact order parameter and 
smooth crossover of deconfinement at high tem¬ 
perature. 

The next remark is about the consistency of 
the quarkyonic regime that was recognized in the 
large-iVc limit and the CSC that is disfavored in 
the large-iVc limit. This is sometimes a source of 


conceptual confusion. To discuss the quarkyonic 
regime we do not have to take the large-fVc limit; 
while the quarkyonic picture was initially stud¬ 
ied in the large-iVc limit, its implication is more 
general and the essence can be discussed even in 
two-color QCD (Brauner et al. 2009; Hands et al. 
2010). If NM and the 2SC phase are identifiable, 
we should consider that this continuity between 
NM and the 2SC phase is nothing but a clear real¬ 
ization of the quarkyonic regime in the real world 
with three colors. This point is emphasized by 
Fukushima (2014) and is precisely the physics pic¬ 
ture that we make full use of in the present work. 

3. Model and parameters 

In the low density region we adopt the EoS 
of nuclear matter, and in this section let us ex¬ 
plain how we compute the EoS of QM at high 
baryon density. We will treat the NJL model 
within the mean-field approximation (for details, 
see e.g. Buballa (2005); Klahn et al. (2007); Kojo 
(2015)). We have three gap energies, A^d, Ads, 
and Asu, associated with diquark condensates. 
Moreover, we should introduce three chiral con¬ 
densates, (uu), (dd), and (ss), which are all dy¬ 
namically determined from the gap equations. As 
we discuss in details, the vector coupling is an es¬ 
sential ingredient for our prescription, so we need 
to treat n = (qj^q) as a mean-field variable, as 
well as three chemical potentials /ig, gs, and /Xg to 
impose electric and color charge neutrality. There¬ 
fore, in total, there are ten mean-field variables 
and we must solve ten gap equations simultane¬ 
ously. 

With these mean-field variables (with a nota¬ 
tion; (/„ = (uu), (pd = (dd), and = (ss)) we can 
write the pressure in the following way: 

~ G's(</u + (/d + ~ 'iGdpu'PdP’s + 

- J^i^l.d + ^ds+^lu) + Pe+P^l , (1) 

where ei{p) represents the energy dispersion rela¬ 
tions obtained from the Dirac Hamiltonian in the 
Nambu-Gorkov doubled basis. The last terms, Pg 
and Pfi, represent the pressure contributions from 
electrons and muons. 
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In the presence of the mean-field n, the quark 
chemical potential is renormalized as = fJ- — 
2G^n. Then, the quark chemical potential takes a 
from of 9 X 9 matrix given by 

(1 = ^J,r^c ® If - Melc ® Qe + (Ms'T'c + ® If > 

( 2 ) 

where Qe = diag(2/3, —1/3, —1/3) is the charge 
matrix in flavor space. The explicit form of the 
gap matrix with Aud, Ads, and Agu can be found 
in the literature (Fukushima ct al. 2005). 

For the model parameters, we use the set of 
Hatsuda & Kunihiro (1994): 

A = 631.4MeV, GgA^ = S.67, GdA^ = -9.29, 
mu,d = 5.5 MeV, m, = 135.7 MeV, (3) 

which are fixed to reproduce hadron phenomenol¬ 
ogy in the vacuum. The other parameters H and 
Gv are left unfixed, and we will use this freedom 
to adjust H and Gv to realize the quark-hadron 
duality interpolating the APR and CFL EoS’s. 

4. Interpolating the EoS with Gv-running 
2SC 

Because the CFL solution does not exist in the 
density region near NM, which indicates a first- 
order phase transition as discussed later, a nat¬ 
ural candidate for the QM ground state there is 
the 2SC state. As elucidated comprehensively in 
Sec. 2, the diquark condensate {ud) in the 2SC 
phase breaks no new symmetry and so it can co¬ 
exist in the hadronic phase. This means that the 
2SC phase could be smoothly connected to NM 
and, furthermore, the quarkyonic scenario requires 
the presence of duality region where NM and 2SC 
can represent the same physical system. 

It is still non-trivial how to formulate such du¬ 
ality in a practical way. We use the APR hadronic 
EoS for the baseline, while the 2SC phase has sev¬ 
eral unconstrained model parameters such as H 
and Gy. Our strategy is to make use of these un¬ 
certainties positively to reproduce APR in the low 
density region in a spirit of the quarkyonic sce¬ 
nario. To this end we treat Gv in the 2SC phase 
as a /iB-dependent control parameter and fit the 
resulting EoS with APR. Besides, we emphasize 
that this choice of Gv as a control parameter would 
provide us with convenient and transparent intu¬ 
ition about modified QM EoS, and at the same 


time, that Gv should be a representative of all un¬ 
known effects like confinement near the normal nu¬ 
clear density. The important point is that in this 
way, once Gv(^b) is determined, we can smoothly 
extend APR toward QM without any artificial ma¬ 
nipulations. 

The pressure P in the 2SC phase depends 
monotonically on Gy, so that we can easily find 
Gv that reproduces APR at each ^b- One might 
think that a /iB-dependent Gy may change one of 
the gap equations, dP/dn = 0. For technical sim¬ 
plicity we do not change the gap equation but treat 
n as an internal mean-field variable which takes a 
different value from riB. Thus, to obtain ub cor¬ 
rectly, we have to compute dP/dp,B including the 
derivative on Gv{fJ-B)- 

In order to find a suitable running coupling, 
we found it necessary to take a diquark cou¬ 
pling larger than 1. Otherwise the NJL pres¬ 
sure at any positive Gy undershoots APR at typ¬ 
ical /is, because the NJL effective quark mass, 
Mq ~ 336 MeV, allows the pressure to develop 
only when > 3Mg ~ 1010 MeV(> rriM) {ruM'- 
nucleon mass). This situation is cured by includ¬ 
ing the sufficiently large diquark coupling which 
effectively reduces the effective quark mass (Kojo 
ct al. 2015). Similar effects can be seen in con¬ 
stituent quark models where the color hyper-fine 
interaction makes the nucleon mass smaller than 
three times constituent quark mass. We will be 
interested in Gy ^ 0.5Gs to achieve sufficient stiff¬ 
ness, and then H > 1.4Gs is required to pass the 
2Mq and causality constraints. Therefore we will 
examine H = 1.5Gs and H = 1.6Gs in the follow¬ 
ing. 

In our fitting procedures one of questions is 
how to choose the fitting range. Let us first con¬ 
sider a relatively wide fitting range, ^b = (1-02 ~ 
1.20) GeV or ub — (1.5 ^ 3.5)no. (We omitted a 
very dilute region /tb ^ 1-02 GeV because it re¬ 
quires too much fine-tuning in our fit functions.) 
We show our numerical results with running Gy 
for H — 1.5Gs (lower solid curve) and H = 1.6Gs 
(upper solid curve) in Fig. 1. We emphasize that 
we did not assume any functional form a priori 
and Gv{p,b ) shown in Fig. 1 results solely from the 
fit to APR once we make a choice of the diquark 
coupling H. We vary H to check the sensitiv¬ 
ity and will see that this choice is near the upper 
limit not to violate the causality. It is important 
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Fig. 1.— Vector coupling in the 2SC phase fitted 
with APR for H = l.SGs (lower solid curve) and 
H = 1.6Gs (upper solid curve). The interpolating 
fit results to the CFL phase with d = 0.4 (See 
Eq.4) are represented by the dotted curves. 



Fig. 2.— Constituent quark masses and the di¬ 
quark gap And in the 2SC phase as functions of 
/tb in the case ot H = 1.5Gv (corresponding to the 
solution shown by a lower solid curve in Fig. 1). 
The chiral phase transition is a smooth crossover 
because of the vector coupling. 


to note that the chiral phase transition is a very 
smooth crossover in the presence of large Gy, so 
that this 2SC phase can accommodate both di¬ 
quark and chiral condensates for any /tb. This is 
clear in Fig. 2 where the constituent quark masses, 
Mu, Md, Ms, and the gap energy A„d are given 
as functions of fiB- There is a small discrepancy 
between Mu and Md because of the electric charge 
neutrality condition that breaks isospin symmetry. 

Our fit to APR tells us that Gy must be very 
large at low and relax to Gy ^ 0.5Gs at large 
^lb- The importance of Gy can be inferred from 
nuclear physics (e.g. w-meson exchange) and typ¬ 
ical magnitude Gy ~ O.SGg in our fit is under¬ 
standable within the conventional context. On 
the other hand, very large value at low rep¬ 
resents something else and requires explanations. 
Let us recall that using the diquark coupling we 
were able to make the onset chemical potential of 
the NJL pressure close to APR, as required from 
the quark-hadron duality. But there remains an¬ 
other problem in our fitting procedures; pressure 
and number density grow much faster in the NJL 
model than in APR. The obvious reason is the 
lack of confinement in our model; a quark excites 
individually, not as a part of a baryon, produc¬ 
ing artificial excess of pressure. But large Gy can 
be used to imitate confining effects by tempering 
the growth of quark number density. This effect 
should be stronger at lower density; as baryons 
overlap, Gy becomes less responsible for confining 
effects and should relax to the value expected from 
the repulsive nuclear forces. 

Interestingly, we have found that the best fit 
form of Gy (/is) is an inverse logarithm for both 
H — 1.5Gs and H = l.GGg. Such an inverse loga¬ 
rithmic is quite suggestive because it is consistent 
with the common form of the running coupling 
constant at one-loop level. However, the valid¬ 
ity of this fitting should be lost at some point of 
the baryon density. In fact, at sufficiently high 
baryon density there is a phase transition from 
the 2SC to CFL phase which softens the EoS. If 
we use the running Gv{^b) determined above, the 
EoS including the CFL phase does not support the 
massive neutron star with M > 2Mq; Gy should 
be ^ 0.5Gs or greater in the CFL phase. Thus 
we need stiffer 2SC pressure with which the EoS 
can remain stiff enough even after the transition 
to the CFL phase occurs. 









The above determination used fitting over the 
range ub — (1-5 ~ 3.5)no, but below we shall relax 
the condition and use only the range ub — (1-5 
2)no, which corresponds to fiq = (340 ~ 345) MeV 
(i.e., fiB = (1-02 - 1.035) GeV). Then the 2SC 
pressure can be stiffer than before in the range of 
ns ~ (2 ~ 3.5)no. It is reasonable to stop fitting 
at fiB = 1.035 GeV, because the baryon density is 
Ub — 2no beyond which APR is no longer reliable. 

To satisfy the boundary conditions, i.e., the 
smooth connection to APR in the lower-density 
side and to the GFL phase with Gv ^ 0.5Gs in 
the higher-density side, we modify Gv(^b) from 
an inverse logarithm to the following form: 


Gv(mb)/Gs 


log[(MB - b)/c] 


(4) 


with an offset by d. Once we fix d, we can de¬ 
termine other three parameters, a, b, c using the 
smooth connection to APR. We changed d to find 
that the massive neutron star with M > 2Mq is 
impossible with d < 0.3. We shall therefore choose 
d = 0.4 throughout this work. The parameters 
fixed in such a way, for H = 1.5Gs and 1.6Gs 
respectively, are listed in Tab. 1 and the corre¬ 
sponding Gv(^b) that interpolates between APR 
and the GFL phase is overlaid by dashed curves in 
Fig. 1. 

With this running-Gv we can find the GFL so¬ 
lution as well as the 2SC phase and then we can 
locate a first-order phase transition between them 
by comparing the pressure. In Fig. 3 we show an 
example for H = 1.5Gv to find a first-order phase 
transition at /ie = 1-31 GeV where the pressure 
of the interpolating 2SC phase with running Gv 
and the GFL phase crosses. We make a remark on 
the connection between APR and the 2SC phase 
around /xb 1 GeV. From Fig. 3 one might think 
that APR has a slightly larger pressure above the 
fitting region, and so APR would be rather fa¬ 
vored. To resolve such confusion we here again 
emphasize our picture of the quarkyonic scenario. 
The change from NM to QM is not any phase tran- 


H/G, 

d 

a 

b [GeV] 

c [GeV] 

1.5 

0.4 

0.05283 

0.4049 

0.5735 

1.6 

0.4 

0.1127 

0.2942 

0.6804 


Table 1: Parameters for the interpolating Gv(/xb) 
between APR and the GFL phase. 


sition but what we assume is a dual regime around 
^b ~ 1 GeV in which NM is gradually taken over 
by QM. In contrast to this smooth crossover from 
NM to QM, the change from the 2SC phase to 
the GFL phase is a genuine physical phase tran¬ 
sition with different symmetry-breaking patterns. 
In many model studies including the present work, 
this phase transition turns out to be of first order. 

Now that we have the EoS for the whole range 
of /xb from NM to CSC, we can compute not only 
the pressure P but the energy e = /xb^b — -P as 
well. Actually, the relation of P vs e is essential for 
the estimation of the neutron star mass. Because 
e involves a first derivative in txb, its value jumps 
discontinuously at the first-order phase transition. 
We can see this behavior in our numerical results 
shown in Fig. 4. It is also clear in Fig. 4 that the 
2SC part hardly changes with different choices of 
H. We can explain this minor dependence from 
the fact that we impose the same boundary condi¬ 
tion of APR at lower density. The other boundary 
condition of the GFL phase side is, on the other 
hand, loosely constrained by APR, and so there 
remains H dependence in the GFL part as is the 
case in Fig. 4. But as we will see below, the value 
of H is constrained by the causality condition. 

One can deduce the sound velocity (squared) 
Cg = dP/de from the slope of curves in Fig. 4. 
We can numerically take the derivative of P with 
respect to e and present obtained in Fig. 5. It 
is important to check the causal condition that 
Cg should not exceed the unity (i.e., the speed of 
light in the natural unit). In the present setup, 
as understood from Fig. 5, the causality is not 
violated. 

The structure of is almost the same irrespec¬ 
tive of the choice of PI ; with increasing e it mono- 
tonically increases in the APR region and then 
has a peak in the 2SG phase. Once a first-order 
phase transition to the GFL phase occurs, is 
pushed down to ^ 0.5. In the limit of large e or 
high baryon density, asymptotically approaches 
^0.6 which slightly depends on P[. An interest¬ 
ing observation in Fig. 5 is that the peak height 
strongly depends on P[ and if P[ is greater than 
^ 1.6Gs, it would go beyond the unity, which 
would violate the causality. Therefore, such a 
large PI is not allowed and H = 1.6Gs is close to 
the upper limit for our choice oid = 0.4 in Eq. (4). 
In this way we find that there is not much uncer- 
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Fig. 3.— Pressure comparison of the 2SC phase 
(dashed curve) and the CFL phase (solid curve) for 
H = 1.5Gs. The pressure of APR (dotted curve) 
is also shown for reference. 



0 0.2 0.4 0.6 0.8 1 1.2 

Energy [GeV/fm^] 

Fig. 4.— Pressure P as a function of the energy 
e for = 1.5Gs (solid curve) and H = 1.6Gs 
(dashed curve). The pressure of APR extrapola¬ 
tion (dotted curve) is shown for reference, though 
the plotted range is outside of its validity region. 


tainty in the EoS determination after all. 

Finally we comment on the connection of our 
result to the pQCD EoS. If we extrapolate our 
P(^_b) curve to 3 GeV, it undershoots the 

pQCD pressure. This is not quite surprising: in 
the pQCD domain, one should expect smaller Gs, 
Gv, and H than those used in our analyses, since 
non-perturbative effects (gluons) are irrelevant by 
definition of the pQCD domain. By simply imple¬ 
menting the reduction of NJL parameters at large 
one can connect our EoS to the pQCD EoS at 
MS > 3 GeV. 

5. Mass-radius relation 

To solve the M-R relation from the TOV equa¬ 
tion what we need is the EoS shown in Fig. 4. 
Plugging our EoS to the TOV equation and chang¬ 
ing the initial condition that is the central pressure 
at r = 0, we can get a curve in the plane of the 
mass and the radius of the neutron star as is just 
the standard procedure. 

We summarize our results in Fig. 6. The M- 
R curves below Mq are essentially determined by 
the APR EoS up to riB 2no. (The tail at 
large R is determined by the crust EoS, for which 
we adopt the SLy model by Douchin & Haensel 
(2001).) Then, the 2SC curves start to take off 
from the APR by degrees. The 2SC curve’s hav¬ 
ing the larger radii than the APR originates from 
a stiffer EoS in the 2SC phase at ub > 2no. In the 
vicinity of M ^ 2Mq the curves have a turning 
point which signals the phase transition from the 
2SC phase to the CFL phase. The maximal mass 
in our parametrization reaches ~ 2.2Mq. 

The neutron star radii at the canonical mass 
1.4M0 typically spread over i? = (9 ~ 16) km, 
depending on the EoS for asymmetric nuclear mat¬ 
ter. The radius of i? ~ 12 km in our results be¬ 
longs to a group of small radii, reflecting that the 
APR at low density is softer than other typical 
hadronic EoS such as mean-field EoS. The small 
radii are consistent with recent observation that 
favors R = (10 ~ 13) km (Ozcl et al. 2010, 2015; 
Steiner et al. 2010, 2013), though those analy¬ 
ses still contain several assumptions to be verified. 
On the other hand there is also a group insisting 
i? > 14 km (Sulcimanov et al. 2011), which fa¬ 
vors stiff hadronic EoS at low density (for hybrid 
EoS with stiff hadronic EoS, see e.g. Bonaimo & 
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Sedrakian (2012); Benic et al. (2015)). Another 
hint for the small radii comes from recent quan¬ 
tum Monte-Carlo calculations which indicate that 
nuclear EoS at low density is relatively soft, or 
even slightly softer than the APR (Gandolfi ct al. 
2012). The heavy-ion data by Danielewicz et al. 
(2002) also favor a soft EoS at low density. 

Our crossover EoS yields relatively simple 
topology for the M-R curves similar to those in 
purely hadronic EoS. In case of the hybrid EoS 
with a first order phase transition, there are more 
varieties of the M-R curves. The classification 
was done by Alford ct al. (2013); Alford & Han 
(2015), assuming the constant speed of sound. 

6. Conclusions 


Fig. 5.— Sound velocity squared as a function 
of the energy for H = 1.5Gs (solid curve) and 
H = 1.6Gs (dashed curve). The APR extrapola¬ 
tion (dotted curve) is shown for reference. 



R [km] 


Fig. 6.— M-R relation for H = 1.5Gs (solid 
curve) and H = l.GGg (dashed curve). The APR 
result (dotted curve) is shown for reference. 


Quarkyonic matter is a likely candidate to 
bridge a gap between nuclear matter and quark 
matter. It is an intermediate state of matter with 
a duality region that is describable by baryons 
or quarks. Thus, the EoS construction is not a 
patchwork of distinct EoS’s of nuclear and quark 
matter, as is a conventional hybrid description 
having first-order phase transitions. 

We discussed the interplay among quarks, di¬ 
quarks, and baryons. Through processes of ex¬ 
changing mesons or quarks between baryons, the 
diquark correlation develops in the intermediate 
state. While the density is low, the diquark is 
combined with other colored excitations to form 
local and color-singlet objects. For this reason 
diquarks should become manifest increasingly as 
baryons involve more and more quark exchange. 
We classified scenarios of deconfinement crossover 
according to the diquark characters and addressed 
what quarkyonic matter implies. 

To interpret the crossover EoS in the language 
of microscopic degrees of freedom, we introduced 
a running vector coupling for quark models, with 
which the quark-hadron duality is implemented in 
a practically feasible manner. The running cou¬ 
pling is constrained by the conventional nuclear 
EoS in the low-density side and by massive neu¬ 
tron stars with M > 2Mq in the high-density side. 
Because the idea is quite generic, one can employ 
any quark models and we adopted an NJL-type 
model for a concrete demonstration in this work. 

The important finding is that the model param¬ 
eters around ne lOng must be comparable with 
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those in the vacuum (i.e. matter at ^ IOuq still 
belongs to a non-perturbative regime); otherwise, 
the 2Mq constraint is not satisfied. The upper 
bound of the diquark interaction comes from the 
causality condition. After all, there is not much 
uncertainty remaining in the parameter determi¬ 
nation. In principle, all the interactions in the 
quark model are mediated by gluons, and so sub¬ 
stantially large couplings require non-perturbative 
gluon dynamics. We need to clarify the origin 
of non-perturbative dynamics in the high-density 
limit in a similar sense to the high-temperature 
limit at which the magnetic sector is still confin¬ 
ing. 

In this paper we interpret the enhancement of 
many-quark (meson) exchanges at ns ^ 2no as 
an indicator for a new state of matter. While we 
consider diquarks in this paper, another branch 
of our picture would be the condensation of 
mesonic objects since many mesons are available 
among baryons. Although the previous studies in 
hadronic models typically predicted strong first 
order phase transitions associated with meson 
condensed phases, recent studies in quark mod¬ 
els indicate that the appearance of inhomoge¬ 
neous mesonic condensates does not necessarily 
accompany strong first order phase transitions, 
thus can avoid the significant softening (Buballa 
& Carignano 2015). The Skyrme crystal descrip¬ 
tions at high density (Pacng et al. 2015) have 
some overlap with such quark descriptions since 
both of them can address several effects beyond 
purely hadronic descriptions, such as the struc¬ 
tural change of hadrons. 

In this work we focused on the bulk properties 
of matter, i.e., the EoS only. To explore the re¬ 
fined characters of matter and to promote the ad¬ 
vantage of our proposed method, we should study 
transport processes such as the cooling problem 
involving the strangeness. We shall leave this for 
future works. 
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No. 15H03652 and 15K13479 (K. F.) and by NSF 
Grants PHY09- 69790 and PHY13-05891 (T. K.). 
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